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ABSTRACT 

We propose a new way of looking at the Baryon Acoustic Oscillations in the Large 
Scale Structure clustering correlation function. We identify a scale Slp that has two 
fundamental features: its position is insensitive to non-linear gravity, redshift space 
distortions, and scale-dependent bias at the 0.5% level; it is geometrical, i.e. indepen¬ 
dent of the power spectrum of the primordial density fluctuation parameters. These 
two properties together make spp, called the “linear point”, an excellent cosmological 
standard ruler. The linear point is also appealing because it is easily identified irrespec¬ 
tively of how non-linearities distort the correlation function. Finally, the correlation 
function amplitude at slp is similarly insensitive to non-linear corrections to within 
a few percent. Hence, exploiting the particular Baryon features in the correlation 
function, we propose three new estimators for growth measurements. A preliminary 
analysis of slp in current data is encouraging. 

Key words: large-scale structure of Universe 


1 INTRODUCTION 


Soon after the revolutionary discovery, through the care¬ 
ful observation of type la supernovae, that the expansion 
of our universe is accelerating (Perlmutter et al.|1999 Riess 


et al. 


several other probes were proposed to verify and 
quantify this important phenomenon. One of the most useful 
employs the so called Baryon Acoustic Oscillations (BAO). 
Acoustic waves propagated through the ionized plasm of the 
early universe, and then were quickly damped at a redshift 
Zd ~ 1050. In linear theory, the baryon distribution retains 
some memory of this process {Peebles & Yujl970). This has 
motivated using the acoustic horizon scale r c j as a cosmolog¬ 
ical standard ruler (e.g. Bassett & Hlozek 20091. The length 
scale is imprinted, for example, in the Large Scale Struc¬ 
ture (LSS) that we measure today ( |Eisenstein et aL||2005[ l; 
it produces a peak in the correlation function (CF) of a va¬ 
riety of tracers of the dark matter distribution at r ~ 100 
Mpc/h (e.g. |Anderson et al.|2014) |Delubac et al.||2015| and 
references therein). 

Like most standard rulers, the acoustic horizon is im¬ 
perfect. Several effects slightly shift or distort the peak in the 
CF, including linear physics (|Seo fe Eisensteinj2dd7 Sanchez| 


* E-mail: stefan o.an sel mi@case .edu 
■ E-maii: gienn.starkmanBcase.edu 

E-mail: shethrkBpliysics.upenn.edu 


et al. 20081, non-linear gravitational evolution (Crocce & 


Scoccimarro 20081, redshift space distortions and scale- 
dependent bias ( |Smith et aL||2008| ). The resulting system¬ 
atic effects are a serious challenge to using the BAO peak 
for sub-percent precision cosmology. Many approaches have 
been suggested to mitigate these effects. (A non-exhaustive 
list includes Crocce^_ScoccimarroJJ2008JJTUiruy^&HiraJ 


matsu|2008||Anselmi fe Pietroni|2012||Desjacques et al. 2010| 
Baldauf et al.||2015| .) 

In this paper we propose a new way of looking at the 
BAO features in the correlation function. We identify a point 
in the correlation function that has three remarkable prop¬ 
erties: 


(1) it is geometric, i.e. its position is independent of the 
power spectrum of the primordial density fluctuation pa¬ 
rameters; 

(2) it is unaffected by nonlinearities and scale-dependent 
bias at percent-level precision, i.e. its position is redshift- 
independent; 

(3) the amplitude of the matter and halo correlation func¬ 
tions at this point are very weakly affected by nonlinearities 
or scale-dependent bias. 

This point lies midway between the scales associated with 
the peak and the dip in the correlation function on BAO 
scales, a definition that holds true even as the dip and the 
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peak shift due to nonlinear physics. Given (2) and (3) we 
denote this the “linear point” (LP). Properties (1) and (2) 
make the LP a standard ruler; property (3) makes the LP a 
useful window into growth of structure and bias. We discuss 
both below. 

We devote most of this paper to an investigation of the 
nonlinear corrections to all the three relevant BAO scales, 
i.e. the dip, the peak and the linear point, taking into ac¬ 
count the dominant effects of nonlinear gravity, redshift- 
space distortions and scale-dependent bias. The three scales 
show different sensitivities to these effects. We find the peak 
to be the most shifted, while the dip is shifted in the opposite 
direction. Due to the near-antisymmetry of the CF across 
the linear point, the position of the LP is much more sta¬ 
ble. This makes the LP a potentially useful standard ruler. 
Furthermore, the matter and halo correlation functions eval¬ 
uated at the LP agree with the linear theory values to a few 
percent. We exploit these features to propose new ways of 
measuring the linear growth and linear bias, and to provide 
consistency checks of our findings with present and future 
data. 

This paper is organized as follows. In Section[2]we intro¬ 
duce the peak, the dip and the “linear point” in the linear 
theory correlation function of dark matter. We show that 
the linear point is almost a purely geometrical point, em¬ 
phasizing how the correlation function is nearly antisym¬ 
metric about this linear point. In Section[3]we introduce the 
main sources of nonlinear distortions that affect the three 
BAO scales, namely nonlinear gravity, redshift distortions 
and scale dependent bias. We quantify systematics from 
nonlinear evolution in Section [4] by showing how the near¬ 
antisymmetry of the CF across the linear point is the key 
feature that drives the linear behavior of the linear point. 
The correlation function amplitude is also linear to a few 
percent at the dip-peak middle point. Section [5] discusses 
a number of possible estimators of the linear point. It also 
leverages the linear theory behaviour of the LP to furnish 
estimates of the linear theory growth factor. A final section 
summarizes our findings and discusses future directions. 


2 THE GEOMETRIC “LINEAR POINT” IN 
THE CORRELATION FUNCTION 

There were two principal ideas behind Baryon Acoustic Os¬ 
cillation (BAO) measurements: 

(i) the comoving baryon acoustic scale at the Compton 
drag epoch corresponds to the baryon acoustic peak s p of 
the matter correlation function; 


(ii) Baryon Acoustic Oscillation measurements are a purely 
geometrical probe, because the baryon acoustic scale con¬ 
tains no information about primordial density fluctuations, 
and so is independent of the initial amplitude A s , the spec¬ 
tral index n s and the optical depth to the last scattering 
surface t. 


It is well known that both these assumptions break down at 
the linear ( |Seo fe Eisensteln||2007| jSanchez et al.||2008[) and 


nonlinear levels (Crocce & Scoccimarro 2008 Smith et al. 


20081 when one probes the BAO at a precision of a few 


percent, which is relevant in this age of precision cosmology. 



Figure 1. BAO scales in the correlation function: the peak, a 
dip and the mid-point between the two, which we call the “linear 
point” (see text). 


While it is well known that s p does not quite equal rd , in 
the following we will identify an actual linear geometric point 
in the correlation function whose position is not affected by 
nonlinear distortions at the percent level: it does not shift 
with time, i.e. redshift. We name this point the “linear point” 
(LP). The linear point is therefore an actual standard ruler, 
unlike the peak scale s p in the correlation function. 

We start our analysis by introducing in Fig. [T] three 
scales in the BAO correlation function: a peak, a dip and 
the “linear point,” defined (for reasons we explain below) as 
the peak-dip middle point: slp = (sd + d p )/ 2. The scale slp 
corresponds to the inflection point in £ to better than 0.5% 
precision; this will be important when we consider scale- 
dependent bias in the next sections. 0 

In Fig. [2] by performing a linear analysis we show that 
the position of the linear point is indeed independent of n s , 
r and A s at the 0.02% level. (The A s independence is true 
by definition in linear theory. Furthermore we checked that 
varying simultaneously n s and r does not change the results 
presented in Fig. [ 2 ]) [j While this geometric property holds 
true in linear theory for the dip and the peak scales as well, 
it is well known that nonlinear evolution affects the posi¬ 
tions of the peak and the dip, spoiling their geometric na¬ 


ture and their redshift independence (Crocce & Scoccimarro 
2008 Smith et al. 2008), and thus their status as standard 


rulers. 


1 In what follows, we focus on an analysis of §. In principle, one 
might have thought that r 2 £ was at least as appropriate. We 
discuss shortly why we have chosen £ instead. 

2 We set the cosmological parameters = 0.0486, Cl c = 0.259 
and Ho = 67.74 to be those of the fiducial model from the |Planck| 
|Collaboration| ( |2015[ ) , and we compute the linear matter power 
spectrum using the CAMB code ( [Lewis et al.|[2Q00| ). While the 
scales Sd, s p and slp depend on these parameters, the ‘geomet¬ 
rical’ nature of spp is generic. 
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Figure 2. Dependence of the “linear point” spp on the primordial density fluctuation parameters, n s , r and A s , expressed relative to 
its value in a fiducial model, which we take to be that of the [Pl anck Collaboration] f2015). Smooth red line shows that spp is independent 
of n s , r and A s to the 0.02% level. 


The linear correlation function is anti-symmetric in the 
peak-to-dip range (see Fig. [I]). One measure of this symme¬ 
try is to compare the scale slp which was defined by the 
peak-dip positions (sd + d p )/ 2, with the scale sa on which 
the amplitude of £ lln equals [£ h71 (sd)+£ h7l ( s p)]/2. These two 
scales differ only slightly: sa ~ 1.002slp- In addition, the 
amplitude £, hn (sLp) corresponds to £ hn (sA) to ~ 2 — 3% 
precision. (We show below that nonlinear effects improve 
this amplitude agreement to 1% level.) This anti-symmetry 
will be crucial for understanding why the linear point slp is 
nearly unaffected by nonlinearities, and is the main reason 
why we prefer to work with £ rather than r 2 £. 

The linear-theory analysis presented in this section thus 
reveals that the matter correlation function gives us a new 
geometric point in the correlation function, the dip-peak 
midpoint, called the “linear point.” In the next sections, we 
show that the linear point is unaffected by nonlinearities at 
percent level precision, making it an excellent cosmological 
“standard ruler.” 


3 NONLINEAR DISTORTIONS 
3.1 Real space 


Nonlinear gravitational clustering smooths the peak in the 
correlation function and shifts it slightly (p3h aradwaj| |1996| 
Seo & Eisenstein 2007 Crocce & Scoccimarro 20081. The 


smoothing is due to displacements from the initial positions, 
so the one-dimensional velocity dispersion in linear theory, 




d 3 g P lm (q,z ) 
(27 r ) 3 q 2 


(1) 


where P hn is the matter power spectrum (PS) in linear the¬ 
ory, sets the scale of the smoothing. 


In particular, in the Zel’dovich approximation (Zel- 
dovich 1970), the smoothing due to the displacements is 
well-approximated by 


P n \k ) 


y -P lln (k) 


( 2 ) 


( |Crocce fe Scoccimarro|[2006| |Peloso et al.||20 j~5] |Vlah et al.| 
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2015). The resulting equation for the nonlinear f is 
d k k 3 P Hn (k) _ fc 2 2 ( 


e\r)«J 


(*) 


k 


2n 2 


jo(kr 


(3) 


Equation has been shown to accurately describe the cor¬ 
relation function in the BAO regime (Crocce & Scoccimarro 
|Angulo et al. | |2008l |Rasera et al.||2014| |Vlah et al. 


2008 


2015 


Peloso et al.|2015 Baldauf et al.|2()l~5 i and it has suc¬ 


cessfully been employed for BAO reconstruction (Eisenstein 
et al. |2007| |White|20l5 1. 


3.2 Redshift space 


Galaxies are observed in redshift space, where the Hubble 
flow induced cosmological redshift has been distorted by pe¬ 
culiar velocities. The real-space-redshift-space map in the 
plane-parallel approximation is 


s = x — fv z z . 


(4) 


Here / = dlnD/dlna is the growth rate, z is the line of 
sight direction, and v(x) = — u/(Hf), where u(x) is the 
peculiar velocity and T-L is the conformal Hubble parameter. 

The same bulk motions which generate nonlinear clus¬ 
tering also give rise to redshift space distortions jK aiser 
19871. Moreover, the angular dependence of equation |4^ 
means that the clustering signal can (and should!) be de¬ 
composed in spherical harmonics. For example, the redshift- 
space monopole associated with £ nl is 


«-<•)■ i/>/ 


d k k 3 P lin (k) 


k 


2n 2 


(1 + M 2 /) 2 


D -fe 2 <r“(l+ M 2 /(2+/)) 


jo(ks ), 


(5) 


where /r = k - z is the cosine between the line of sight and the 
wave vector. The integral over fi can be performed analyt¬ 
ically, making numerical integration over k very fast (Des- 
jacques fe Sheth|2010[ Peloso et al.|2015 1. The linear theory 
limit has no nonlinear smearing; so, setting a v —> 0 yield s 
£o(s) = (l + 2//3 + / 2 /5) £(s) in linear theory (Kaiser] 1987). 
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3.3 Halo bias: linear and scale dependent 

In the real Universe (as opposed to simulations thereof), we 
cannot directly observe the matter field. Instead, we measure 
galaxies that inhabit halos which, in turn, are biased tracers 
of the matter field. In this section we focus on how halo bias 
(the halo-matter relation) affects the correlation function in 
the BAO regime. A detailed analysis of how bias affects the 
BAO features of the correlation function is beyond the scope 
of this paper. However, we present a preliminary investiga¬ 
tion of this topic, and expect to carry on this study in future 
works. 


We use the peaks theory approach to halo bias (Bardeen 


ct al. 19861. As before, we employ a simplified description 


that allows us to focus on the main effects (e.g., we ignore 
mode coupling terms). A good approximation to the corre¬ 
lation function of peaks is 


,nl,, fdkk 3 P Hn (k) _ fe 2 CT 2 
x [bio(z) + boi(«)fc 2 ] jo(kr), 


( 6 ) 


where bf 0 is the usual Eulerian linear bias (Desjacques et al. 
|2010 1. Similarly, the redshift-space monopole is simply 


£S,nl , y 

M),/i h\ a > 2 




/ 


d k k 3 P lin {k) 


k 


27r 2 


[&io i z ) + b 0 


(z)k 2 + /r 2 /] 


s -fe 2 ^(i+M 2 /(2+/)) j- 0 ( fcs ) ( 7 ) 


where the integral over fi can again be done analytically 
( Desjacques fe Sheth}201 0). 

Strictly speaking, we should replace a v in the expres¬ 
sions above with the velocity-biased value appropriate for 


peaks: a p < a v (Desjacques & Sheth|20l0). In addition, the 


expressions above should include the smoothing windows of 
radius R p (one for each peak) associated with identifying 
the peaks. However, these are almost degenerate with a p in 
real space (and not quite as denegerate in redshift space). In 
effect, the expressions above have set R p + a 2 ~ a 2 . In this 
respect, we compute a v from equation 0 while we treat 
equations 0 and 0 as effective formulae, where bf 0 , bo\ 
could be evaluated by fitting to simulations (as was done in 
|Peloso et al.||2015 l. 


4 WHY DOES THE LINEAR POINT REMAIN 
LINEAR? 

4.1 Dark matter 


Nonlinearities affect the locations of the BAO peak, dip and 
LP scales in two different ways: 


(a) smoothing the correlation function; 

(b) shifting the correlation function 

(e.g. |Smith et al.| [2008|. To see this clearly, rewrite equa¬ 
tion 0 as 

rV)=r‘(ixi;ifc) (8) 


where Rg = v2cr„ and 


r'(|x|;R) = |dV^_^l^ ( | x '| ) . (9) 



Figure 3. Monopole of the redshift-space correlation function 
rescaled by the linear Kaiser factor of (1 -\- 2//3 + / 2 /5). Vertical 
lines show the scales of the peak, dip and linear point for linear 
theo ry (d ashed), the ID Gaussian convolution (dotted) (i.e. equa¬ 
tion |lOl with r'/r — > 1), and the nonlinear model £q"( r) (contin¬ 
uous) . 


The 3D convolution can be simplified to 


G(|x| ;R)= J dr' 


(r-r') 
/ „- 7777T 


(r+r ') 2 


r e 
r 


2 R 2 — e 2 R 2 


(r - r ') 

~ I dr' 1 - 6 2R2 


/■ 


r ( 2 ttR 2 ) 1 /2 


(27rR 2 ) 1 / 2 

j-lin / /\ 

i ( ?- ) ■ 


j-lin / /\ 

C (r) 


( 10 ) 


where the last equality is practically exact for r > 2 R. This 
shows explicitly that the exponential term inside the integral 
smooths the correlation function, while the r'/r term shifts 
the peak of the Gaussian filter to higher r. This, in turn, 
moves the peak of £ to smaller r. 

Similarly, equation 0 for the redshift-space monopole 
can be written as 


£o’T) = ^ f^{i +l x 2 f) 2 e\\^SG) (ii) 


where Sg = \/l + /r 2 /( 2 + /). Since Sg > Rg the 

smoothing is larger in redshift space, but the overall effect 
is otherwise similar to that in real space. 

To illustrate the effect of these nonlinearities, Fig. [3] 
shows the monopole of the redshift-space correlation func¬ 
tion rescaled by the linear theory factor of D 2 (1 + 2//3 + 
f 2 / 5) at z = 0 in a ACDM cosmology with = 0.265, 
£h,h 2 = 0.0226, h = 0.71, n = 0.963 and a 8 = 0.80 The 
vertical lines show the positions of the peak, dip and linear 
point for linear theory (a v —> 0) (dashed), the ID Gaussian 
convolution (i.e. equation (101 with r'/r —> 1) (dotted), and 


3 Note that these are not the Planck cosmological parameters we 
used in Section [2] rather, they are those of the simulation set we 
use later in the paper (see Section [40. This is why, the peak, dip 
and LP scales are ~ o//~ 1 Mpc larger than in Fig. Ill 
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the nonlinear model £q l (r) (continuous). Notice that the ID 
Gaussian convolution moves the peak and the dip by the 
same amount, but in opposite directions, whereas the r'/r 
term shifts the whole correlation function to smaller values 
of r by 1%. Thus, the two effects add for the peak but they 
subtract for the dip. As a result, the dip shifts less than the 
peak. 

More importantly, the scale of the linear point shifts 
even less than either the peak or the dip. The net result is 
that the peak moves by 3.5% to the left, the dip by 1.5% to 
the right and the linear point by just 1% to the left. (These 
shifts are smaller at higher redshifts, of course.) The peak 
is consequently the most sensitive point to nonlinear correc¬ 
tions and the linear point is the least. Moreover, because the 
dip and the peak share the same nonlinear smoothing, the 
amplitude of the linear point remains unchanged to within a 
few percent even though the dip and the peak (heights and 
scales) change significantly. 


4.2 Biased tracers: Halos 


To understand how halo bias affects the correlation function 
in the BAO regime we focus, for simplicity, on the real space 
results. The same arguments apply to the redshift space 
monopole. 

To gain insight we write equation |6| as a convolution 
of the linear correlation function, and simplify as before to 
obtain 


CUr) * 


/ 


_ (r— t / ) 2 

r ' p 2 R% 

dr' — £ lln ( r ') -_ 

r 5 1 1 {2-kR 2 g Y/i 


x nl 0 + 2b w boi Rb (r r) 


R% 

+ &§i (r-r') 4 - 6( r r'^-l- 3 ^ ,,^) 


where we have dropped the 
venience. 


R% 


‘E” on the bias factors for con- 


Desjacques et al. (20101 have emphasized the fact that 


the terms proportional to b oi and &oi contribute significantly 
to the shape of £hh on BAO scales. As their Figure 3 shows, 
these terms are approximately anti-symmetric in r, so they 
serve to enhance the BAO feature by increasing the height 
of the peak and depressing the dip. However, their figure 
also shows that these terms contribute little to the scale or 
the amplitude of the correlation function at the linear point 
- a point they do not discuss. 

There are two ways to see why this happens. First, note 
that the term oc ben scales as k 2 in Fourier space, so it is 
cx V 2 £ !irl = (2 /r)d£ lin /dr + d 2 £ lin /dr 2 in real space. On 
the peak and dip scales, the first term is zero, but the second 
term contributes maximally to the amplitude of ^hh- How¬ 
ever, at the inflection point, the second term is zero by defi¬ 
nition, and the first term is almost surely very small on BAO 
scales because biased tracers assemble their mass from re¬ 
gions which are much smaller than the BAO scale. Hence, at 
the inflection point, the net contribution from the boi term is 
negligible. The argument for the 6^ term is similar: It scales 
as k 4 in Fourier space, so it is oc V 2 (V 2 £) = 4£®/r + 

(i.e., to third and fourth derivatives of £ ilrl wrt r ). Since 
£ ! * n is rather smooth at the inflection point, these higher 


derivatives are small, and so the net contribution from the 


term 6 qi term is small. Figure 3 in Desjacques et al. (2010) 
shows that, indeed, while the terms oc &oi and fejji con¬ 
tribute maximally around the BAO peak and dip scales, 
they are substantially smaller at the inflection point. Since 
slp lies very close to the inflection point, we expect the 
scale-dependent bias terms to be strongly suppressed at slp, 
making ££ l h (s L p ) « b \ 0 ("'(sifjQ 

Alternatively, we can modify the convolution argument 
of the previous subsection to see how the terms related to 
&oi and 6 qi contribute to £(^. To begin, we note that 


/ 


dr' 


(r-r'F 


(2nR%y/ 2 


— (r — r') 2 

~R l 


= o 


(13) 


and 


(r — r ') 2 


dr' 


( 27 tr 2 g : 


1/2 


(r - r') 4 - 6(r - r') 2 R 2 G + 3 R 4 C 
R% 


= 0 


(14) 

also. If £ hn (r) is approximately antisymmetric in r, then 
the integrals above will remain very small even if we insert a 
factor of £ ln (r') in the integrands. The terms proportional 
to foot and boi include an additional r'/r factor, but this 
factor will not change the value of the integral dramatically. 
We have checked numerically that including it indeed makes 
little difference. As a result, the terms proportional to &oi 
and &o i contribute little to the result, in agreement with our 
discussion in the preceding paragraph (and with Figure 3 of 
Desjacques et al!|2010l. Thus, this provides another way to 
see why £^ ( slp ) is insensitive to scale-dependent bias even 
when feoi is not negligible, despite the fact that the height 
and scale of the peak and dip can be strongly modified. 

We conclude that the antisymmetric structure of the 
correlation function is preserved even in the presence of 
scale-dependent bias. While this scale-dependence moves the 
peak and the dip in opposite directions, the linear point re¬ 
mains fixed. Furthermore the amplitude of the biased cor¬ 
relation function at the linear point is only sensitive to the 
scale-independent part of the linear bias factor, bio- Both 
these reasons increase the appeal of the LP as a standard 
ruler candidate. 


4.3 Comparison with simulations 

To check the accuracy of our analysis, we first compare equa¬ 
tions § and to measurements in the N-body simulations 
of a ACDM cosmology with fi m = 0.265, fibh 2 = 0.0226, 
h = 0.71, n = 0.963 and as = 0.8 carried out by |Sato~^z| 
Matsubara ( 2011|). 


The top panel of Fig.[4]shows the real-space correlation 
function at four different redshifts, rescaled by the scale- 
independent linear growth factor. The bottom panel shows 
the same rescaled quantity for the redshift-space monopole. 


4 There is another inflection point just beyond the BAO peak- 
scale. Therefore, we expect the scale-dependent bias terms to be 
suppressed there as well, and this is indeed seen in Figure 3 of 
|Desjacques et ah] ( |2010| l . While exploiting this larger scale as an 
additional standard ruler is beyond the scope of this work, we 
return briefly to this scale in the next section. 
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s [h 1 Mpc] 


Figure 4. Top: Real space correlation function rescaled by the 
linear growth factor at four redshifts. Bottom: Similarly rescaled 
redshift-space monopole. Smooth lines show equations ([3]) and 
; symbols with error bars show the N-body results. The vertical 
dashed red line marks the scale of the linear point in linear theory. 


The symbols with error bars show the measurements while 
the continuous lines show the theoretical predictions pre¬ 
sented in the previous section. (We compute these as inte¬ 
grals over P(k) rather than as real-space convolutions.) 

The simulation measurements scatter much less than 
the 5 — 10% one would expect given the error bars, which 
were obtained from the scatter among 30 different realiza¬ 
tions. Apparently, the points are strongly correlated, and a 
full error analysis must properly treat the covariance. If the 
smooth curve that the eye would fit through the points is a 
good approximation of the correlation function to be inferred 
from the simulations, and the very small point-to-point scat¬ 
ter around that by-eye fit is a fairer approximation of the 
true error bars, then it appears that, while the amplitude of 
the correlation function is at times discrepant between the 
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Figure 5. Top: Real space BAO peak, dip and LP scales. Bot¬ 
tom: Redshift space BAO peak, dip and LP positions. Smooth 
black curves show the measurements in simulations, shaded gray 
regions show 1% variation around the measurements, and dotted 
and dashed curves show the linear and nonlinear theory values of 
these scales. 


simulations and the theory, the scales of the dip and the peak 
are correctly predicted to subpercent accuracy, and the LP 
to 0.5%. Fig. [5] shows this explicitly: the solid curves show 
the BAO peak, dip and LP scales, which we estimate by 
fitting a high-order polynomial to the measured £, and the 
corresponding dotted and dashed curves show the linear and 
nonlinear theory values of these scales. 

Fig .[6] shows a similar analysis of the clustering of halos 
more massive than > 1.37 x 10 12 /i _1 A/q at all the red- 
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Figure 6. Same as Fig. [4] but for biased tracers (halos more 
massive than > 1.37 X 10 12 h _1 A7g). 
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shifts. [^] Symbols with error bars show the N-body results 
and smooth curves show eqns. 0 and 0- For these, the lin¬ 
ear bias bi Q (z) is fit to the simulation, and we found &oi(z) 
to be negligible given the large errorbars. The top panel 
shows the real-space correlation function at four different 
redshifts, rescaled by the scale-independent linear growth 
factor as well as the linear bias bj 0 (z). The linear bias val¬ 
ues at z = (0,0.5,1,2) are b 10 ~ (1,1.29,1.72,2.95). The 
bottom panel, shows the rescaled redshift-space monopole, 
where /3 = f /b. The nonlinear theory predictions are in rea¬ 
sonable agreement with the measurements. Again, the error 
bars appear to be highly correlated, and, while the ampli- 

5 |Sato &: Matsubara| |201l[ | identified the halos in each output 
using a friends-of-friends (FOF) group finder with linking length 
set to 0.2 times the mean particle separation; the mass limit above 
corresponds to halos with more than 20 particles. 


Figure 7. Same as Fig. [5] but for biased tracers. 


tude of £ is often discrepant by up to 10%, the locations of 
the standard rulers are much more robust. Fig. [7J the ana¬ 
log of Fig. [5] shows a more quantitative comparison: Notice 
that, especially in redshift space, the peak and dip scales dif¬ 
fer more significantly from the measurements. Nevertheless, 
the LP scale still agrees to better than 0.5%. We conclude 
that the LP is indeed a good ruler even for biased tracers at 
late times in redshift space. 

Before we exploit this fact, it is worth making one fi¬ 
nal point. Recall that the amplitude of £ ltn (sLP ) is within 
about 3% of [£ h "(sd) + £ iln ( s p)]/2- Figures [ 4 ] and [6] show 
that nonlinear evolution makes £ nl less steep than £ hn , both 
in real and redshift space. Hence, £ nl (sLp) corresponds to 
[{ n («d) + 5™*(sp)]/2 at percent level precision, a factor of 
two improvement with respect to linear theory. 
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5 TWO CONSEQUENCES 

5.1 Towards more accurate distance 
measurements 


We have made the point that the LP is mainly sensitive to 
the shift arising from the r'/r factor (see Section 4.1l. We 
have found that a simple redshift independent 0.5% correc- 
tion removes this effect: 


slp 


Sp Sd 


x 1.005 


(15) 


and restores agreement with the linear prediction to 0.5% at 
all redshifts both in real and redshift space. 

This is not the only option. Since the peak and the 
dip scales are geometrical in linear theory, any linear com¬ 
bination of Sd and s p which cancels the nonlinear shifts will 
return a geometric LP as well. We have found that the op¬ 
timum weighted sum at z = 0 is 


sun,w — 0.7s d + 0.3s p ; 


(16) 


this combination returns the same linear theory value at 
all higher redshifts to 0.7% (and can be corrected to 0.35% 


using the same approach as in equation (151. 


One could even go further and employ a redshift 
weighted sum 


sun(z) = c(z)s d + d(z)s p 


(17) 


where c(z) and d(z) could be extracted from simulations so 
that sun(z') would correspond to the linear theory predic¬ 
tion, i.e. 1/2 < c(z ) < 7/10 and 3/10 < d(z) < 1/2. But we 
have not pursued this further as we feel we have made our 
main point: using information from both the peak and the 
dip is potentially very powerful. Indeed not only do eqns. 


(151 (171 provide very stable scales for distance measure¬ 


ments, but comparing them at different redshifts provides 
useful cross checks. 


5.2 Using the LP to measure growth 

So far we have made the point that the LP is a good standard 
rod. However, its amplitude is also rather well-understood, 
so we now explore its use in measuring the linear theory 
growth factor. 

Over the BAO range of scales, our model predicts that 
the peak and dip heights should differ more from their linear 
theory values than do the peak and dip scales (c.f. Figures [4] 
and[6|. However, at the LP, the amplitude of the correlation 
function is rather well described by linear theory. To make 
this point, Fig. [8] shows the ratio of the nonlinear to the 
linear CF in the simulations (symbols with error bars) and 
in our model (smooth curves). Although our model is able 
to describe the scale dependence of this ratio rather well 
(changes from unity can be as large as 40% at 2 = 0), it is 
remarkable how little the amplitude of the LP changes with 
redshift. Fig. [9] shows a similar analysis of biased tracers. 
Again, at the LP, this ratio is within a few percent of unity]^] 


6 As we anticipated in footnote 4 at the end of Section 4.2, there 
appears to be a scale larger than the BAO peak scale on which 
this ratio is close to unity at all redshifts even for biased tracers. 
While it is not as close to unity as is the ratio on scale sjp, this 



r [h 1 Mpc] 


2.0 



s [h 1 Mpc] 


Figure 8. Ratio of the nonlinear correlation function to that 
in linear theory in real (top) and the redshift (bottom) space at 
four different redshifts. Smooth curves show the analytic model, 
symbols with error bars show the measurements and the vertical 
dashed red line shows the scale of the LP in linear theory. At the 
LP, this ratio is very close to unity at all times. 


For unbiased tracers, our analysis of the LP scale mo¬ 
tivates the following estimator for the ratio of the linear 
growth factor at redshifts z and z': 

P 2 G) _ Z hn (s L p,z) £{slp,z) , . 

D 2 (z’) e in (SLP,Z') £(s' lp ,z>) ' 

Here slp = ( s p + Sd)/ 2, where carated quantities are the 
estimates inferred from observations. 

While this is conceptually straightforward, it assumes 
errors on the measured amplitude are negligible. Since this 


scale and the amplitude of § on it potentially allow additional 
consistency checks of our spp methodology which we leave for 
future work. 
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Figure 9. Same as Figure [8] but for biased tracers (halos more 
massive than M d > 1.37 X IQ 12 h 1 A/q ). 


estimators of the ratio 

D\z) 1 + f f{z) + \f\z) 


D 2 (z')l + lf(z') + lP(z') 

from measurements of the redshift-space monopole. 

Extending this to biased tracers is similarly straightfor¬ 
ward. We have already argued that the scale dependent bias 
terms matter little, so inserting the reds h ift s pace monopole 
of biased tracers for £ in equations (181 (201 furnishes esti¬ 
mates of 


b 2 10 (z)P 2 (z) l+lP{z)+\0 2 {z) 
b 2 o{z')D 2 {z') 1 + 1/3(2') + j8 2 (z') 


( 21 ) 


(The corresponding sums over the real space correlation 
function yield estimates of this quantity with /3 —> 0.) 

In all three versions of our estimator, the shift of the 
whole correlation function due to the r'/r term in the con¬ 
volution (see above) is irrelevant, provided that we consis¬ 
tently compare the peak-dip middle point[^] Therefore, these 
three estimators could also be used as a consistency check 
for the antisymmetry of the correlation function. In ACDM, 
we expect this antisymmetry to be slightly broken by the 
mode-coupling contribution to the correlation function, al¬ 
though exploring this is beyond the scope of this paper. 

While implementing the three estimators with actual 
data is beyond the scope of this paper, we believe that the 
following example is encouraging. We focus on the halo red- 
shift space monopole at 2 = 0 and z' = 0.5, since these 
biased tracers at low redshift are where the differences from 
linear theory are expected to be largest. We then simply 
compare the theoretical predictions of OAhh. with measure¬ 
ments in simulations (obtained by fitting a polyno¬ 

mial to the measured £) and with the linear theory value. 
For the halos in the simulation, the linear theory value of 
the ratio defined in (211 is 0.9347. Our nonlinear model pre¬ 
dicts 0.9359 for all three estimators of this ratio, whereas 
the measured values are (0.9303, 0.9338, 0.9321). The agree¬ 
ment between the estimators as well as the nonlinear model 
with the linear theory value is better than 1%. Clearly, an 
accurate analysis with actual BAO data is in order. 


is rarely the case, one could reduce the effect of measure¬ 
ment errors by averaging the measured amplitude over more 
scales. E.g., one might set 


P 2 (z) e in (s P ,z)+e in (s d ,z) _ i(8 P ,z) + j(a d ,z) 

D 2 (z') t lin {s p ,z')+£ lin (s d ,z') £{%,,*)+ £(s' d ,z')' 

(19) 

or even 


D\z) 

D 2 {z') 


f s / d dxi lin {x,z ) 
d* £ lin (x, z 1 ) 


E £{xi,z)/N(z) 

Sd< x i<Sp 

E i(xi,z')/N(z')’ 
( 20 ) 


where N is the number of measurements done within the 
peak and the dip. 

In view of our discussion of the similarity of nonlin¬ 
earities in real and redshift space, one could build similar 


6 DISCUSSION AND CONCLUSIONS 


Since the first detection of the Baryon Acoustic peak in the 
galaxy correlation function ( Eisenstein et al.|2005 l, tremen¬ 
dous effort, both theoretical and observational, has been de¬ 
voted to reducing systematic uncertainties and improving 
the statistics. 

With this in mind, we have searched for a special point 
in the correlation function whose scale and amplitude are 
least affected by nonlinear evolution. We have found that 
the point which lies midway between the peak and the dip 
of the correlation function on BAO scales, which we call the 
linear point slp (see, e.g., Figure[lJ, is shifted in scale from 
its linear theory value by less than 0.5%. Previous work has 
shown that nonlinear evolution smears out the height and 
shifts the scale of the peak. We showed that there are sim¬ 
ilar effects on the dip, but that these nonlinear corrections 


We ignore the correction introduced in equation (11 51). 
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cancel out almost exactly at the linear point (Figure [3]). We 
provided analytic arguments for why this same LP is special 
for unbiased as well as biased tracers, in real and redshift 
space (Section 4), and showed that this was indeed the case 
in simulations (Figures [4] and [d|. 


We believe that this motivates use of the LP for dis¬ 
tance measurements, as it may remove the need to fit broad¬ 
ened correlation function templates to the measurement in 
order to infer the acoustic scale as is currently done (e.g., 
Sanchez et al.| [2009 |San chez et al. 2012]). For example, es¬ 


timating slp from the published BOSS DR11 galaxy cor¬ 
relation (Figure 10 of [An derson et ah] [i2014|, following the 
same simple procedures we used to estimate it in the simu¬ 
lations presented in the main text (fitting a high-order poly¬ 
nomial to the measured £) and applying equation ( |15| ) yields 
95.1/i _1 Mpc. For comparison, slp in the fiducial model of 
Anderson et al. ( 2014| ), divided by a = 1.03 as their Ta¬ 
ble 7 suggests, also yields 95.1/i -1 Mpc. We leave a more 
careful implementation of this procedure (e.g., accounting 
for covariances between bins, reconstruction, etc.) to future 
work. 


A nontrivial consequence of the cancellation of nonlin¬ 
ear effects is that it is not just the scale of the LP which 
remains fixed: its height is also given by linear theory, to 
percent level precision, even at late times, for unbiased and 
biased tracers in real and redshift space (Figures [8] and [9]. 
We exploit this interesting property to motivate three esti¬ 
mators of the linear theory growth equations ( |18| |-( [20| ) - 
and to compare our findings with data. 

In practice, our simplest estimators of the growth con¬ 
strain the combination bio D rather than bio and D sep¬ 
arately. We expect the width of the ‘nonlinear’ smoothing 
to be approximately independent of bio, so it may be that 
combining our LP-derived constraint on bio D with the more 
traditional estimators of the acoustic scale and its broaden¬ 
ing will allow one to infer bio and D separately. Breaking 
this degeneracy is the subject of work in progress. 


Our detailed results for the CF structure in the BAO 
regime could be modified by non-standard physics. At the 
linear level we expect that the presence of massive neutrinos, 
some form of dark energy, or some modification of grav¬ 
ity will change the discussion of Section [5] ( |Thepsuriya fe| 


Lewis 2015). On the other hand, when we deal with non- 


linearities we should be careful to make some distinctions. 
The BAO dip-peak region in the presence of massive neutri¬ 


nos is well described by the same physics as ACDM (Peloso 


et al. 20151. Quintessence and clustering quintessence also 


show the same nonlinear propagator as the standard case 
(Anselmi et al. 2014), the transition from large to small 
scales being smooth. On the other hand, a preferred scale 
often appears in modified gravity models, so equation ©> 
can no longer be applied naively (Scoccimarro 2009). It is 
not obvious that the antisymmetric features of the corre¬ 
lation function will be preserved in general. In this regard, 
precise measurements by future surveys of departures from 
the antisymmetry in the BAO range of scales would be a 
clear signature of non-standard gravity. 
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